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ABSTRACT 

The  basic  operating  principles  of  an  experimental  system  for  the 
direct  and  continuous  measurement  of  solid  propellant  burning  rates  are 
presented.  The  system  has  been  designed  and  fabricated.  Several 
components  of  the  measurement  system  have  been  modified  in  order  to 
increase  the  precision  of  the  burning  rate  measurements. 

A  continuous  burning  rate  measurement  technique,  termed  the  Servo¬ 
mechanism  Technique,  will  be  employed  for  obtaining  erosive  burning  rate 
data  for  types  BDI  and  BUU  double-base  propellants.  Burning  rate 
measurements  will  be  made  with  different  gas  flow  velocities  parallel 
to  the  burning  propellant  surface;  the  gas  velocities  will  approach 
the  acoustic  speed. 

A  feasibility  study  was  conducted  for  determining  the  adaptability 
of  microwave  techniques  to  the  measurement  of  the  burning  rate  of  a 
solid  propellant.  It  was  concluded,  because  of  the  dependence  of  the 
microwave  attenuation  upon  the  combustion  conditions  present  in  a 
research  rocket  motor,  that  microwave  techniques  are  not  readily  adaptable 
to  such  burning  rate  measurements. 

A  feasibility  study  indicates  that  a  technique  employing  ultra¬ 
sonic  pulses  can  be  developed  for  obtaining  direct  measurements  of 
the  burning  rate  of  a  solid  propellant.  The  technique  is  based  on 
measuring  the  time  for  an  ultrasonic  pulse  to  travel  through  a  pro¬ 


pellant  sample. 


INTRODUCTION 


Accurate  knowledge  of  the  burning  rate  of  a  solid  propellant  is 
necessary  for  reducing  the  costly  full  scale  testing  of  newly  designed 
solid  propellant  rocket  motors.  To  date,  burning  rate  theories  cannot 
predict  with  sufficient  accuracy  the  burning  rate  characteristics  of 
solid  propellants.  Consequently,  burning  rate  data  which  are  employed 
for  designing  solid  propellant  rocket  engines  must  be  obtained  experi¬ 
mentally. 

The  conventional  experimental  techniques  employed  currently  for 

determining  the  burning  rate  of  a  solid  propellant  often  yield  erroneous 

results  (l)^  Several  of  the  indirect  measurement  techniques  utilize 

some  form  of  probe  which  is  embedded  in  a  propellant  sample  of  propellant 

grain.  Such  techniques  generally  produce  inaccurate  burning  rate  data 

because  the  high  thermal  conductivity  of  the  probe  increases  the  rate  of 

heat  transfer  from  the  combustion  zone  to  the  propellant  surface. 

Furthermore,  embedded  probes  grossly  affect  the  combustion  gas  flow 

field  adjacent  to  the  burning  propellant  surface  which  is  undesirable 

2 

in  stu  ies  concerned  with  measuring  erosive  burning  rates.  Direct 

1.  Numbers  in  parentheses  indicate  references  in  the  Bibliography. 

2.  Erosive  burning  is  defined  as  the  change  in  the  linear  burning  rate  due 
to  combustion  gas  velocity  parallel  to  the  burning  surface.  Linear 
burning  rate  is  defined  as  the  rate  at  which  the  burning  surface  recedes 
when  there  is  no  combustion  gas  velocity  parallel  to  the  burning  surface. 
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methods  of  burning  rate  measurement,  such  as  the  photographing  of  the 
receding  burning  surface  (2),  yield  reliable  burning  rate  data  and  do 
not  disturb  the  flov  field  but  the  reduction  of  the  data  is  tedious  and 
time  consuming.  Moreover,  available  transparent  vindov  materials  char 
quickly  under  high  temperature,  high  gas  velocity  flov  conditions. 

High  temperature  gas  flows  are  unavoidable  in  erosive  burning  rate 
experiments,  especially  in  the  case  of  high  energy  solid  propellants. 

Investigations  at  the  Jet  Propulsion  Center,  Purdue  University, 
have  been  concerned  with  improving  the  techniques  for  measuring  the 
burning  rates  of  solid  rocket  propellants  (2), (3,4,5)?  The  subject 
report  presents  the  recent  progress  in  the  aforementioned  investigations. 

Section  II  of  the  subject  report  is  concerned  with  a  burning  rate 
measurement  technique  which  involves  a  servomechanism,  and  also  dis¬ 
cusses  a  research  program  based  on  employing  the  servomechanism  technique 
for  obtaining  erosive  burning  rate  data  of  solid  propellants. 

Section  III  of  this  report  discusses  the  results  of  a  feasibility 
study  for  adapting  microwave  techniques  to  burning  rate  measurements. 

Section  IV  discusses  a  technique  for  applying  ultrasonic  waves  to 
the  continuous  measurement  of  the  burning  rate  of  a  solid  propellant* 


. .  . ,  .... — i _ _ — . . .  . 

3.  Funded  by  the  Air  Force  Office  of  Scientific  Research  under 
Grant  AF-AFOSR  2O7-63. 


II  SERVOMECHANISM  MEASUREMENT  TECHNIQUE 

An  experimental  system  for  the  direct  and  continuous  measurement 
of  solid  propellant  burning  rates  has  been  developed  at  the  Jet  Pro¬ 
pulsion  Center.  A  detailed  discussion  of  that  technique  Is  presented 
In  References  3,  4,  and  5*  Only  a  brief  description  of  that  servo¬ 
mechanism  technique  Is  presented  here. 

Figure  1  Illustrates  schematically  the  essential  features  of  the 
servomechanism  measurement  technique  for  obtaining  the  burning  rate  of 
a  solid  propellant.  A  sample  of  a  solid  propellant  Is  bonded  to  a 
rectangular  holder  vhlch  is  attached  to  a  shaft  so  that  It  can  move  the 
propellant  sample  In  a  linear  direction;  the  holder  Is  incorporated 
Into  a  t vo- dimensional  research  rocket  motor.  As  the  solid  propellant 
burns,  the  burning  surface  recedes  In  the  direction  of  Its  normal 
with  the  normal  velocity  r  (6).  Simultaneously,  the  servomechanism 
moves  the  propellant  sample  In  the  direction  opposite  to  that  of  the 

receding  burning  surface,  with  the  velocity  r  .  Thus,  if  the  shaft 

8 

velocity  r  equals  the  burning  rate,  r,  the  surface  of  the  propellant 
sample  Is  maintained  at  a  fixed  position  in  the  two- dimensional  rocket 
motor. 

In  that  case,  the  direct  measurement  of  the  velocity  of  the  shaft, 

r  ,  yields  the  Instantaneous  burning  rate  of  the  propellant  sample. 
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PROPELLANT  POSITIONING 


FIG.  I  BURNING  RATE  MEASUREMENT  BY  THE 
SERVOMECHANISM  TECHNIQUE 
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Feedback  System 

The  major  problem  encountered  in  developing  the  aforementioned 
servomechanism  technique  vas  concerned  with  developing  a  surface 
position  detector  vhlch  provides  an  accurate  feedback  signal  to  the 
servomechanism.  Figure  2  presents  one  successful  type  of  feedback 
system f the  Gamma  Ray  Feedback  System.  Another  method  for  providing 
a  feedback  signal  to  the  servomechanism  is  discussed  in  Section  III 
of  this  report. 

The  Gamma  Ray  Feedback  System  (see  Fig.  2)  employs  a  collimated 
beam  of  gamma  rays  vhlch  is  attenuated  by  the  propellant  sample  located 
inside  the  tvo-dimens ional  rocket  motor.  The  gamma  rays  are  provided 
by  a  125  millicurle  source  of  Cobalt *60.  As  the  propellant  bums,  the 
intensity  of  the  collimated  beam  increases  due  to  the  difference  In 
the  density  of  the  solid  propellant  and  its  combustion  products.  A 
scintillation  probe  measures  the  change  in  the  Intensity  of  the  beam 
and  provides  a  feedback  signal  to  the  servomechanism.  The  servo* 
mechanism  then  moves  the  propellant  sample  to  maintain  the  burning 
surface  at  a  predetermined  position. 

To  protect  personnel  from  radiation,  the  radioactive  isotope  source 
is  surrounded  by  lead  shielding.  In  addition,  the  scintillation  probe 
is  provided  with  a  lead  shield  for  reducing  the  background  radiation 
count;  a  high  background  radioactive  count  vould  produce  a  feedback 
signal  that  is  not  related  to  the  propellant  position. 


PROPELLANT  POSITION NG 


FIG.  2  GAMMA  RAY  FEEDBACK  SYSTEM 
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Servomechanism  System 

Figure  3  presents  a  schematic  diagram  of  the  integrated  servo¬ 
mechanism  system.  The  intensity  of  the  collimated  beam  of  gamma  rays 
is  monitored  by  the  scintillation  probe  which  produces  a  pulsed  output. 
The  pulses  of  the  scintillation  probe  are  counted  by  the  ratemeter  which 
converts  the  count  rate  to  a  millivolt  signal.  The  millivolt  output  of 
the  ratemeter  is  amplified  and  subtracted  from  the  reference  voltage, 
and  the  resulting  error  signal  is  amplified  by  the  servo  amplifier. 

The  amplified  error  signal  is  for  exciting  the  field  of  the  rotating 
amplifier;  the  latter  is  driven  at  constant  speed  and  produces  the 
power  for  operating  the  D.  C.  servomotor,  which  operates  with  a  constant 
field  voltage,  and  the  servomotor  repositions  the  propellant  sample 
within  the  rocket  motor. 

The  burning  rate  of  the  propellant  is  monitored  by  two  separate 
transducer  systems.  A  linear  potentiometer,  connected  to  the  shaft 
which  positions  the  propellant  sample  within  the  motor,  produces  a 
signal  which  is  proportional  to  the  position  of  the  shaft;  the  latter 
is  recorded  on  a  direct  writing  oscillograph,  and  the  slope  of  the 
position-time  trace  is  equal  to  the  burning  rate  of  the  propellant 
sample. 

The  servomotor  tachometer  coupled  with  the  oscillograph  provides 
a  system  for  monitoring  the  servomotor  speed.  Thus,  by  knowing  the  gear 
ratio  of  the  propellant  feed  mechanism,  the  burning  rate  of  the  propel¬ 
lant  sample  may  be  read  directly  from  the  oscillogram. 


POSITION 
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The  chamber  pressure  in  the  rocket  motor  is  recorded  simultaneously 
with  the  aforementioned  transducer  signals.  Consequently,  at  each 
instant  during  an  experiment  the  burning  rate  of  the  propellant  sample 
can  be  correlated  with  the  chamber  pressure. 

Improvement  of  Servosystem 
Propellant  Feed  Mechanism 

As  reported  in  Reference  5;  the  method  for  moving  the  propellant 
sample  vithin  the  rocket  motor  vns  found  to  be  inadequate  under  certain 
operating  conditions.  At  chamber  pressures  above  400  psi  the  counter 
torque  vhich  the  propellant  holder  exerted  on  the  propellant  position¬ 
ing  shaft  was  found  to  be  excessive.  Consequently,  the  propellant 
holder  had  a  tendency  to  Jam  between  the  sides  of  the  rocket  motor 
causing  the  servomechanism  to  become  inoperative.  The  propellant  feed 
mechanism  was  redesigned  to  overcome  the  aforementioned  difficulty,  and 
to  reduce  the  friction  of  the  propellant  positioning  mechanism. 

Figure  4  illustrates  schematically  the  redesigned  propellant  feed 
mechanism.  To  provide  a  torque  to  counter  the  friction  torque  produced 
by  the  internally  threaded  worm  gear,  a  "T"  bar  was  attached  to  one  end 
of  the  propellant  positioning  shaft.  The  "t"  is  supported  at  each  end 
by  parallel  oil  impregnated  bronze  bars.  Thus,  as  the  propellant  posi¬ 
tioning  shaft  moves  either  in  a  forward  or  reverse  direction  the  sliding 
"T"  bar  provides  the  necessary  counter  torque  for  preventing  rotation  of 


the  shaft. 


WORM  GEAR — i 

TEST - 1  |  r - BRONZE 

ROCKET  \  LINEAR - »  A  PARALLEL 
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BALL  THRUST  BEARINGS 
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Figure  k  also  shows  the  redesigned  bearing  system  incorporated 
into  the  propellant  feed  mechanism.  To  reduce  shaft  friction,  linear 
hall  hearings  (Thomson  Industries)  have  been  installed  instead  of 
bronze  journal  hearings.  Furthermore,  the  bronze  thrust  hearings 
employed  formerly  were  replaced  hy  ball  thrust  hearings . 

The  redesigned  shaft  feed  mechanism  has  been  tested  at  chamber 
pressures  in  excess  of  1,000  psi,  and  the  propellant  feed  mechanism 
performed  satisfactorily  under  those  conditions. 

Isotope  Source  Strength 

The  statistical  error  (root  mean  square  error/count  rate)  encountered 
in  monitoring  the  Intensity  of  the  collimated  gamma  ray  beam  is  a  function 
of  the  counting  rate  and  the  time  constant  of  the  ratemeter.  The  net 
result  of  the  statistical  error  is  a  fluctuation  in  the  feedback  signal 
to  the  servomechanism;  the  latter  can  be  reduced  to  a  minimum  (excluding 
electronic  noise)  by  either  of  two  means:  1.  The  time  constant  of  the 
ratemeter  may  be  increased,  or  2.  Tie  nominal  counting  rate  at  the  desired 
propellant  position  may  be  increased  by  increasing  the  activity  of  the 
isotope  source. 

Increasing  the  time  constant  of  the  ratemeter  to  reduce  the  fluctua¬ 
tions  of  the  feedback  signal  to  the  servomechanism  results  in  an  overall 
decrease  in  the  performance  of  the  integrated  system.  For  the  subject 
servomechanism,  optimum  time  constant  settings  have  been  found  to  be 
0.1  second  or  less  (k).  To  reduce  the  feedback  signal  fluctuations  to 
an  acceptable  value  (employing  a  50  millicurie  Cs-137  isotope  source) 
time  constants  greater  than  j.O  seconds  are  necessary.  Thus,  since  a 
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J.O  second  tine  constcnt  ic  detrimental  to  overall  system  performance 
the  er.l\  alternative  to  reduce  the  feedback  signal  fluctuation  vas  to 
increase  the  activity  of  the  isotope  source. 

The  law  governing  the  feedback  fluctuations  can  be  expressed,  for 
all  practical  purposes,  by  the  so-called  Poisson  error  distribution  (7). 
Employing  that  concept,  the  feedback  signal  statistical  error  encountered 
in  men iter in  ;  th  *  collimated  Lean  intensity  is  given  by  the  following 


cruet ion.  Thus, 


E-  C 
n 


i/a 


n 


(1) 


where  E  is  the  feedback  signal  statistical  error,  <T*  is  the  root  mean 
square  error,  and  Ti  is  the  nominal  counting  rate  cs  indicated  by  the 
ratemeter . 


Ac  reported  tn  R.  fcrcnce  the  nominal  counting  rate  when  5°  nc 

cf  C.  -I77  wr.c  employed  as  a  source  of  gamma  radiation  was  85,000  counts 
per  minute.  However,  at  a  ratemeter  time  constant  setting  of  0.1  second, 
85,000  counts  per  minute  is  equivalent  to  approximately  142  counts  per 
0.1  second.  At  that  nominal  counting  rate,  an  elementary  calculation 
shows  the  statistical  error  to  be  in  excess  of  8$. 

To  reduce  the  statistical  error  the  intensity  of  the  collimated 
beam  was  increased  to  approximately  2.75  x  10^  counts  per  minute,  by 
replacing  the  50  nc  Co-1  7  isotope  source  by  125  me  of  Cobalt-60.  Thus, 
at  the  optimum  ratemeter  time  constant  settings  the  statistical  error 
is  found  to  be  lesn  than  2$. 
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To  monitor  the  increased  intaasity  of  the  collimated  gamma  ray  beam 
it  was  necessary  to  replace  the  original  ratemeter  (Baird- Atomic  Model  432-A) 
by  one  having  a  higher  counting  rate  capability  (Baird-Atomic  Model  Cs-400). 
With  the  modified  isotope  source  the  servomechanism  system  has  operated 
satisfactorily. 

4 

Erosive  Burning  Rate  Study 

The  servomechanism  technique  for  measuring  the  burning  rate  of  a 
solid  propellant  is  being  applied  to  the  investigation  of  erosive  burning 
characteristics  of  both  double-base  and  composite  solid  propellants.  A 
previous  erosive  burning  study  at  the  Jet  Propulsion  Center  employed  a 
photographic  technique  as  the  means  of  obtaining  burning  rate  data  (8). 

The  principal  disadvantage  of  the  aforementioned  photographic 
measurement  technique  was  that  the  reduction  of  the  photographic  data  is 
tedious  and  time  consuming.  Also,  available  transparent  window  materials, 
with  which  the  research  rocket  motor  was  equipped,  chared  quickly  under 
high  temperature, high  gas  velocity  flow  conditions.  The  application  of 
the  servomechanism  technique  for  erosive  burning  rate  measurement  elim¬ 
inates  the  aforementioned  disadvantages  because  the  burning  rate  is 
recorded  directly  by  an  oscillograph,  and  no  windows  are  needed  in  the 
test  section  of  the  research  rocket  motor. 

Figure  5  illustrates  schematically  the  experimental  apparatus  to  be 

employed  for  obtaining  erosive  burning  rate  data.  A  gas  generator  will 

ITT  This  work  will  be  co- funded  by  AFOSR  under  Grant  AF-APOSR  207-64  and 
by  Ball  is  tic  Research  Laboratories,  Aberdeen  Proving  Grounds,  under 
Contract  DA-ll-022-AMC-2088(x) . 
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FIG.  5  EROSIVE  BURNING  RATE  APPARATUS 
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be  used  to  supply  hot  combustion  gases  to  a  tvo-d linens ional  test  section; 
the  latter  will  contain  a  propellant  sample  which  will  be  moved  by  the 
servomechanism  at  a  rate  equal  to  its  burning  rate.  Since  one  side  of 
the  two-dimensional  flow  channel  in  the  test  section  will  be  slightly 
divergent,  the  Mach  number  at  any  cross-sectional  reference  plane  in 
the  test  section  will  remain  constant  during  the  experiment. 

The  initial  erosive  burning  rate  studies  will  be  conducted  with 
types  BDI  and  BUU^  double-base  propellants.  Future  studies  will  use 
double-base  and  composite  propellants  with  additives  such  as  ballistic 
modifiers  and  different  percentages  of  aluminum;  in  that  manner,  the 
effect  of  the  additives  upon  the  erosive  burning  rate  will  be  determined. 

The  Mach  number  of  the  combustion  gases  in  the  reference  plane  of 
the  test  section  will  be  maintained  at  0,  0.2,  0.4,  0.6,  0.8,  and  0.95 

with  test  section  static  pressures  of  600,  800,  1000,  and  1200  psi. 

Tie  initial  temperature  of  the  gas  generator  propellant  and  the  propel¬ 
lant  sample  will  be  -40,  70,  and  140°F.  It  is  worth  noting  that  in 
the  past  it  has  been  difficult  to  make  erosive  burning  rate  measurements 
at  Mach  numbers  above  0.6.  The  servomechanism  measurement  technique 
facilitates  making  measurements  at  high  Mach  numbers. 


5.  Complete  propellant  specifications  may  be  found  in  the  SPIA/M2 
Propellant  Manual,  the  Johns  Hopkins  University,  Silver  Spring, 
Maryland,  Confidential. 
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Progress  of  the  Erosive  Burning  Rate  Investigation 
Teat  Facilities 

To  provide  adequate  facilities  for  performing  the  erosive  burning 
rate  experiments  a  tect  cell  at  the  Jet  Propjision  Center  has  been  equipped 
with  the  Instrumentation  for  operating  the  servomechanism  system  and  instru¬ 
mentation  for  measuring  aM  recording  pressures.  Included  in  the  instrumenta¬ 
tion  are  means  for  the  remote  operation  of  high-speed  Fastax  cameras. 

Gas  Generator  Grain 

The  design  of  the  gas  generating  propellant  grain  has  been  completed. 
Propellant  grains  and  propellant  samples  for  the  subject  study  will  be 
supplied  (GFE)  by  the  Ballistic  Research  Laboratories,  Aberdeen  Proving 
Ground,  Aberdeen,  Maryland. 

The  grain  geometry  which  was  chosen  for  the  gas  generator  is  a  slotted 
tube  configuration  with  two  axial  slots  at  the  aft  end  of  the  grain  (9)- 
The  slotted  tube  configuration  which  was  chosen  has  the  distinct  advantage 
of  an  essentially  constant  raasc  flow  rate  (burning  area  change  of  less 
than  2%)  while  requiring  a  simple  mandrel  configuration. 

Igniter 

A  modified  basket  type  igniter  (Aerojet-General,  Model  CASI-Il)  will 

be  used  to  ignite  the  gas  generating  propellant  grains.  The  main  charge 

of  the  igniter  will  consist  of  approximately  50  grams  of  ALCK)  pyrotechnic 

pellets  and  the  charge  in  it  ia  ter  will  be  a  cartridge  type  squib  (Holex 

Model  1 395> )  *  Free  volume  tests  of  the  igniter  are  planned  in  order  to 

predict  approximate  igniter  performance. 

£7  Trade  Name  of  Aerojet-General  Corporation. 
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III  BURNING  RATE  DETERMINATION  BY  A 
MICROWAVE  TECHNIQUE 

This  section  discusses  the  results  of  a  feasibility  study  for 
determining  the  adaptability  of  microwave  techniques  to  the  measurement 
of  the  burning  rate  of  a  solid  propellant.  It  should  be  mentioned  at 
the  outset,  that  currently  the  complex  processes  Involved  In  the  attenua¬ 
tion  of  microwaves  by  solid  propellant  combustion  products  are  relatively 
unknown.  The  analysis  which  follows  is  at  best  a  first  order  approximation 
to  the  complex  phenomenon  of  microwave  attenuation  in  solid  propellant 
combustion  gases. 

Figure  6  illustrates  schematically  a  system  for  providing  a  feedback 
signal  to  the  subject  servomechanism.  Its  operating  principle  involves 
passing  a  microwave  beam  through  the  test  rocket  motor  and  measuring  the 
change  in  the  effective  intensity  of  the  received  beam  as  a  function  of 
the  position  of  the  movable  propellant  sample. 

The  small  windows  in  tht  sides  of  the  chamber  are  made  of  a  di¬ 
electric  material,  such  as  quartz  or  Plexiglas,  which  provide  a  high 
transmittance  fcr  the  microwaves.  To  be  capable  of  predicting  the 
characteristics  of  the  proposed  feedback  system  it  is  necessary  to  examine 
the  propagation  characteristics  of  microwaves  in  the  combustion  products 
of  the  solid  propellant  sample. 


h* 
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The  Microwave  Equation 

A  mlcrovave  is  a  type  of  electronagnetic  wave  which  propagates  at 
the  speed  of  light  vith  frequencies  between  1000-300,000  megacycles  per 
second.  Its  wavelength,  however,  is  longer  than  that  for  either  ultra¬ 
violet,  visible,  or  infrared  radiation.  Microwaves  are  emitted  and 
absorbed  during  the  transitions  within  the  molecular  rotational  spectra; 
i.e.,  the  lov  molecular  energy  levels  (10). 

The  propagation  of  the  microwave  can  be  expressed  by  the  following 

7 

general  form  of  the  wave  equation.  Thus 


V  (<7  •  E)  -  E 


(2) 


where 


E  =  the  electric  field  intensity,  volt/m 
i  =  the  convective  current,  coulomb/ m2  sec 
|i  =  the  y  meability  of  the  medium,  henry/m 

€  =  the  permissivity  of  the  medium,  farad/m  *  coulomb  /joule  m 
v  =  the  gradient  operator,  + 


Further, 

V-  7  = 


(3) 


where  is  the  charge  density  in  coulomb/ft\  To  solve  equation  2, 
it  is  assumed  that  the  wave  is  not  subjected  to  attenuation.  In  that 

p 

case  p  and  i  both  vanish  from  equation  2.  Equation  2  reduces  to 
6 


V2  E  =  M€ 


(4) 


7.  A  derivation  of  this  general  wave  equation  is  presented  in  Appendix  A. 

8.  Refer  to  equation  30  in  Appendix  B.  Since  pg  ■  0;  hence,  ^  *  E  »  0. 
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The  solution  of  equation  k  for  a  one-dimens ionel  w nvr  propagation  is  of 
the  form 

E  A  -  ,<Z)  .  B  c‘!(uH  *  K7-'  (5) 


where 


A,B  are  integration  constants  to  be  determined  by  the  appropriate 
boundary  conditions 
K  is  defined  as  K  wfir,  l/ cm 
^  is  the  microwave  rnuian  frequency,  l/sec 
t  is  the  time,  sec 


is  the  thickness  of  the  combustion  chamber,  cm 


is  defined  as  J  -  1 


If  the  microwave  is  propagated  one  dimensionally  along  the  positive 
Z  coordinate  only,  the  second  term  in  equation  5  can  be  omitted. 

Furthermore,  at  time  equal  to  zero,  and  Z  equal  to  zero  (the  window- 
gas  interface)  the  electric  field  intensity  is  Eq.  Considering  one¬ 
dimensional  wave  propagation,  equation  5  becomes 

E  -  E  '  KZ)  (6) 

o 

The  attenuation  of  the  microwave  is  governed  by  Beer's  Lav  which 
states  that  the  attenuation  of  the  microwave  is  an  exponential  function 
of  the  propagation  distance  Z.  Thus, 

e  .  eo  ^  *  E)  (r) 

where  ot  is  an  absorption  coefficient  tfaich  is  a  function  of  the  medium  in 


which  the  microwave  propagates. 
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The  Absorption  Coefficient 

The  absorption  coefficient  at  is  a  function  of  several  modes  of 
absorption;  conductive,  molecular,  and  electronic-transition,  and  other 
effects  such  as  ionized  combustion  gases.  The  effects  of  ionized  com¬ 
bustion  gases  are  discussed  in  a  Jeter  section  of  this  report. 

The  conductive  absorption,  o^,  may  be  found  by  assuming  that  the 
combustion  gases  are  dielectric,  and  further  for  simplicity  that  the 
charge  density  is  uniform.  According  to  Ohm's-  Law 

T  =  V*  E  (8) 

where  O'"  is  the  electrical  conductivity  in  mho/meter.  Further,  con¬ 
sidering  the  condu  L i  t'or^t’or  n.l  .  emotion  7  can  be  rewritten  in 
the  form 


E  =  E  e 
o 


,Qiz  -  KZ) 


(9) 


where  0^  is  the  conductive  absorption  coefficient  which  is  a  part  of  the 
total  absorption  coefficient  or.  Employing  equation  8,  the  first  term  of 
equation  2  may  be  rewritten  in  the  form 

V  (  V  •  e)  V  [v  •  (-^)] 

’rl'7  v' T)1  U°> 

Combining  equation  ^  and  10,  one  obtains 


V(V  'E) 

*'5-[  ^(VCe>] 


=  0 
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since  it  has  been  assumed  that  the  charge  density  is  uniform  throughout 
the  propagation  medium. 

Substituting  equation  8  into  equation  2,  and  noticing  that  the  first 

9 

term  of  equation  2  vanishes,  the  general  wave  equation  becomes 

2 

-7  2  E  =  -  ix  *£  (<T  E)  -  p€  £-| 

dt 

or 

V?  E  =  U  ♦  ne  1  (11) 

Further,  utilizing  equation  9  and  taking  its  Laplacian,  first  derivative, 
and  second  derivative  with  respect  to  time,  respectively,  one  obtains  the 
following  three  expressions. 


&  = 
at 


2 

E  =  =  (  -  CL  -  JK)2  E 

a  z 

(j  e 


E 


Substituting  these  expressions  into  equation  11  yields 

-  (-  0^  -  JK)?  E  =  -JuTuT'E  +  E 

or 

(c^  ♦  JK)P  =  J  uT  u> 

=  UVT-  K2 


9*  The  vector  symbols  are  deleted  here  since  we  have  assumed  one-dimensional 

propagation  of  the  microwave. 
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hence, 


o^  +  jk  =  k2  +  j  ur  LitrJ 


1/2 


“t  *  - 


] 


1/2 


JK 


JK 


[ 

[‘ 


!  .  J 


1/2 


M* 


- 

ta,*j 


1/2 


For  a  small  loss  due  to  dielectric  attenuation;  i.e.,  1,  the 

above  relation  can  be  expanded  into  a  Taylor  series.  Considering  only 
the  first  terra  then 

1  <T  ^ 

♦  JK  =  JK  (l  -  ^7) 

rK 


JK  ♦ 


2&re 


hence, 


„  -  K 

°!l  2ure 


CT~  U>ln€ 
2  UJ  € 


(12) 


Thus,  the  conductive  absorption  coefficient,  a^,  is  seen  to  be  a 
function  of  the  electrical  conductivity,  permeability,  and  permissivity 
of  the  propagation  medJum. 

To  illustrate  the  magnitude  of  the  conductive  absorption  coefficient, 
consider  the  conductive  absorption  coefficient  for  air  at  2500°K  and 
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500  psi:  air  is  considered  here  because  of  the  lack  of  data  for  the 
electrical  conductivity  of  solid  propellant  combustion  gases.  The 
electrical  conductivity  of  air  at  these  conditions  is  (ll) 

<T"  *  10  4  mho/cm 

The  permissivity  of  air,  is  (12 ) 

€  ”  3S  *  10  ^  farad/m 
The  permeability  of  air,  i3  (12) 


U  =  4n  *  10  henry/m 

Employing  equation  12  to  calculate  the  conductive  absorption  coefficient, 
one  obtains 


0.006 


or, 


=  0.0106k 

=  0.027 


-1 

cm 


-I 

cm 

inch 


Thus,  since  the  attenuation  of  the  microwaves  in  the  combustion  gases 
isafunctioc  of  e  0(2  it  may  be  concluded  that  the  attenuation  of  microwaves 
by  conductive  absorption  is  not  the  dominant  mechanism  of  attenuation. 

To  determine  the  molecular  absorption,  Q^,  that  phenomenon  must  be 
studied  from  a  microscopic  point  of  view. 

Consider  a  microwave  which  is  incident  on  a  combustion  gas  molecule. 
When  the  microwave  collides  with  the  gas  molecule,  an  electron  is  excited 
from  a  quantum  state  1  to  another  quantum  state  ra.  Employing  the  concept 
of  molecular  excitation,  the  molecular  absorption  coefficient  has  been 
shown  to  be  (10) 
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8n  \  N 
3OOO  he 


^lm  no  I  ^lm | 


vhere 


(13) 


N  =  Avogadro's  Number, 

h  =  PlaLch's  Constant, 

c  =  the  velocity  of  light, 

f.,  =  the  mean  frequency  of  quantum  states  1  and  m, 

_Lm 

n  =  the  molecular  concentration, 

o 

subscript  i  refers  to  the  i  molecule  in  the  combustion  gases, 
and  ImJ  is  expressed  by  the  integral, 


(14) 


where 

L  ,  L  are  the  molecular  wave  functions  at  states  1  and  m 
1  m 

respectively,  and 

p  =  the  dipole  moment  of  the  molecule . 

The  evaluation  of  equation  14  is  based  on  the  assumption  of  the 

transition  of  a  partis,  in  a  box  which  makes  possible  the  determination 

of  the  allowable  energies  and  the  position  probability  function  for  a 

particle  moving  one-dimens ionally  in  a  region  of  length  a.  Moreover,  the 

wave  function  is  equal  to  zero  at  x=0  and  x=a.  Thus,  a  is  the  lowest 

common  multiple  of  l/f,  and  l/f  .  The  equation  for  p,  can  be  further 

x  m  im 

simplified  to  (10) 

ea  cos(l  -  m)  rt  -  1  _  cob(1  ♦  m)  n  -  1 

^lm  2  /.  \2  t .  \2 

n  (1  -  m)  (1  +  ra) 


(15) 
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where  e  is  the  electronic  change.  The  molecular  absorption  follows  the 

Rules  of  Section;  that  is,  the  transition  between  energy  levels  can  occur 

only  if  one  level  has  an  odd  quantum  number  and  the  other  an  even  quantum 

number,  as  is  shown  by  equation  15*  For  example,  when  1  =  1,  m  -  3, 

equation  15  yields  u.  =0  and  hence  cr_  -  0.  Since  the  combustion  products 

are  composed  of  several  gaseous  species,  the  molecular  absorption  for  each 

species  is  evaluated  by  applying  equation  15*  The  summation  of  the  a  's 

>  i 

is  the  molecular  absorption  coefficient  for  the  combustion  products.  Thus, 
the  molecular  absorption  coefficient  is  given  by 
n 

a2  ’  Z  “2)1  <l6> 

i  =  l 

where 

<%,  =  the  molecular  absorption  coefficient  for  the  combustion  products, 

a0  =  the  molecular  absorption  coefficient  for  the  i*'*1  species  of  the 

combustion  gas,  and 

n  -  the  total  number  of  the  molecular  species  in  the  combustion 
products . 

The  molecular  absorption  coefficient  is  a  function  of  the  combustion 

gas  temperature  and  pressure.  For  different  combustion  temperatures  there 

exist  different  quantum  states  for  transition;  hence,  f.  and  u  have 

lm  lm 

different  values.  Furthermore,  for  different  combustion  gas  temperatures 
and  pressures  the  molecular  concentrations,  nQ,  is  a  variable.  Consequently, 
the  molecular  absorption  coefficients  are  functions  :>f  the  propellent 
composition  as  well  as  the  combustion  environment. 
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It  is  necessary,  therefore,  to  select  an  optimum  microwave  frequency 
for  each  propellant  composition  and  combustion  pressure.  The  optimum 
frequency  can  be  determined  by  using  different  microvave  frequencies 
through  experimental  trials,  or  by  analyzing  the  spectrum  of  the  radiation 
of  the  propellant  flame  which  will  enable  one  to  select  those  particular 
frequencies  which  are  subject  to  a  minimum  attenuation. 

The  absorption  due  to  electronic  transition  is  usually  within  the 
region  of  high  frequencies;  i.e.,  in  the  visible  or  ultraviolet  regions  (10). 
Thus,  since  the  frequencies  of  the  microwave  are  lower,  this  mode  of 
absorption  can  be  neglected. 

From  the  foregoing  analysis  it  is  concluded  that  the  absorption 
coefficient,  a  (see  equation  j)  is  given  by 

a  =  +  a2 

or 

n 

amal*  Z  *>,1  W 

i  =  l 

Although  the  form  of  equation  17  is  relatively  simple,  the  evaluation 
of  the  absorption  coefficient  Of  encounters  formidable  difficulties. 

Equation  17  requires  that  the  composition  of  the  combustion  gases  near 
the  burning  surface  be  known.  To  date,  the  composition  of  the  combustion 
above  the  surface  of  burning  solid  propellants  cannot  be  predicted  with 
an  accuracy  required  for  determining  a_  . 
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Influence  of  the  Ionized  Combustion  Gases 
The  combustion  gases  in  the  test  rocket  motor  will,  to  a  certain 
degree,  be  ionized  ( 13 ) .  Hence,  the  microwave  will  be  propagated  through 
the  slightly  ionized  gases  which  will  have  the  characteristics  of  a  plasma. 
It  is  to  be  expected,  therefore,  that  the  plasma  will  block  the  propaga¬ 
tion  of  the  microwave  at  certain  frequencies. 

The  attenuation  of  the  microwaves  by  the  plasma  can  be  predicted  by 
applying  Newton's  Law  of  motion.  Thus 


dV 
ra  d^ 


e  E 


(18) 


where  e  is  the  electronic  charge,  V  is  the  velocity  of  the  electron,  and 
m  is  the  mass  of  the  electron. 

Defining  p  as  the  charge  density  and  assuming  the  microwave  is 
e 

propagated  only  in  the  positive  z  direction,  then 


i 


V 


or 

V  =  —  -  —  E  (19) 

pe  pe 

In  other  words,  the  velocity  of  the  electron,  V,  is  also  subject  to  an 
exponential  attenuation  since  it  is  proportional  to  E  at  a  certain 
temperature  and  pressure;  hence, 


f  -  J  °rv 

j  u^v  =  -  E 

0  m 


or, 


).n 
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E 

ml 


thus 


Iml  °e 


1  '  ‘  E 


(20) 


By  assuming  uniform  p  ,  the  wave  equation,  equation  2,  can  be  reduced  to 

-  (a  ♦  JK)?  E  =  -  p  p  E  ♦  m€  e 

m  e 


(a  ♦  JK )‘ 


2  lei 

-  OJ*  ♦  P  -1—-  p 
m  e 


2  2 
=  -  ur  p€  ♦  p€ 


where. 


■t 


(21) 


Equation  21  defines  the  so-called  radian  plasma  frequency,  l/sec. 
Furthermore , 


(a  +  JK)?  -  -  ur2  pe  (1  -  ^|) 

U sd 


Equation  22  shows  that 


a  ♦  JK  =  real 

or,  K  =  0 

In  that  case,  the  microwave  will  not  pass  through  the  ionized  gas. 


Ca/?*  uX ,  a  +  JK  =  imaginary 

a  =  0 
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In  the  latter  case,  the  microwave  will  pass  through  the  ionized  gas  but, 
its  phase  angle  will  be  altered. 

The  degree  of  ionization,  x  =  n^/n^ ,  can  be  calculated  by  the 
application  of  the  law  of  mass  action  and  Nc-rnst's  hoot  theorem  (14). 
Thus 


2 

x 


P 


2.4  x  10 


-eV/kT 

e 


(23) 


where 


n^  =  the  ionized  molecular  concentration  of  the  gas, 

n  =  the  molecular  concentration  of  the  gos, 
o 

p  =  the  partial  pressure  of  a  gas  in  mm  Mg, 

o 

T  =  the  temperature  of  the  gas  in  K, 
eV  =  the  electronic  potential,  and 
k  =  the  Bolt  man  constant. 

Further,  the  chi  rge  density  can  be  calculated  from 


P 


e 


(24) 


where 

Po  =  the  density  of  the  gas, 

M  =  the  molecular  weight  of  the  gas, 

N  =  the  Avogadro's  Number,  and 
e  --  the  cl  arge  of  the  electron. 

As  an  example,  the  plasma  frequencies  for  the  combustion  gases  of 
an  Aerojet-General  Corp.  propellant,  ALT  l6l,  is  shown  in  Table  1.  Th< 
necessary  data  were  taken  from  Reference  6  and  the  electronic  potentia 
from  Reference  14.  All  dsta  presented  in  Table  1  are  based  on  a  mean 
combustion  pressure  of  500  psia  and  a  propellant  flame  temperature  of  1 


Gas 

Mole 

Fraction 

Partial 

Pressure 
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Table  1 

Electronic 
Voltage  (ev) 

Degree  of 
Ionization 

X 

Plasma 

Frequency 
f  ( Gc/sec) 

b2 

.274 

137-0 

15.  b 

I.69  x 

lo'10 

2-55 

h20 

.099 

49-5 

12.6 

2.19  X 

lo'8 

I8.3 

CVI 

o 

o 

O.O33 

lo .  5 

13-7 

6.95  X 

io"9 

1.2 

N2 

0.001 

0.5 

15-5 

2.54  X 

10"  9 

15.8  Mc/sec 

CO 

0.440 

220 

KCJ 

0.151 

75-5 

S02 

0.002 

1.0 

An  important  consideration  In  adapting  a  microwave  detection  system 
as  a  servo  feedback  transducer  is  that  the  plasma  frequency  is  a  function 
of  the  mole  fraction,  partial  pressures,  and  temperature  of  the  combustion 
gas.  Hence,  as  the  burning  rate  varies,  so  does  the  plasms  frequency  and 
accordingly  the  microwave  attenuation.  Consequently,  a  separate  microwave 
transducer  calibration  system  would  be  required  for  adapting  the  trans¬ 
ducer  to  the  servomechanism  measurement  systerr  for  each  different  propel¬ 
lant  and  combustion  pressure. 
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Conclusions 

Based  on  the  foregoing  simplified  analysis  it  is  concluded  that  the 
microwave  detection  system  is  not  readily  adaptable  as  a  feedbacK  trans¬ 
ducer  for  the  servomechanism  measurement  technique,  because  of  the  strong 
dependence  of  the  microwave  attenuation  upon  the  combustion  conditions  in 
the  research  rocket  motor. 

In  view  of  the  above,  no  further  work  is  planned  involving  the 
application  of  the  microwave  technique  to  the  measurement  of  the  burning 
rate  of  a  solid  propellant. 
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IV  BURNING  RATE  DETERMINATION  BY  AN  ULTRASONIC  TECHNIQUE 


This  section  discusses  the  use  of  an  ultrasonic  technique  for 
determining  the  burning  rate  of  a  solid  propellant  sample  which  is  located 
within  a  research  rocket  motor.  The  ultrasonic  method  under  investiga¬ 
tion,  termed  the  pulse-echo  technique,  shows  promise  as  a  technique  for 
the  continuous  measurement  of  the  burning  rate  of  a  solid  propellant. 
Before  discussing  the  application  of  the  pulse-echo  technique  to  burn¬ 
ing  rate  measurements  a  brief  discussion  will  be  presented  of  certain 
characteristics  of  ultrasonic  wave  propagation. 


Propagation  Characteristics  of  Ultrasonic  Waves 
Ultrasonics  is  a  type  of  sonic  pressure  wave  propagated  at  frequencies 
above  the  audio  frequencies  of  approximately  20  kilocycles  per  second. 

The  ultrasonic  wave  is  propagated  at  the  speed  of  sound  but  with  shorter 
wavelengths.  The  intensity  ratio  of  the  reflected  and  incident  waves  at 
a  normal  angle  of  incidence  at  a  boundary  between  two  media  can  be  cal¬ 
culated  by  (15) 


I_  [Z1  ~  Z2 

h  LZ1+Z2 


(25) 


where 


I, I  are  the  intensities  of  the  reflected  and  incident  waves 
respectively,  and 

^,Z0  are  the  specific  acoustic  impedenc®  of  media  1  and  2  respectively. 


The  specific  acoustic  impedence  can  be  calculated  by  (15) 

Z  =  pV  (26) 

where 

p  =  the  density  of  the  medium,  and 

V  =  the  sonic  velocity  of  the  wave  in  the  medium. 

Table  2  presents  the  specific  acoustic  impedence  of  several  media  (16). 


Table  2 


Specific  Acoustic  Impedence  of  Several  Material 


Material 


Specific  Acoustic  Impedence  x  10-  e] cm  sec 


Equation  25  can  be  employed  for  calculating  the  percentage  of 
reflected  energy  between  two  dissimilar  materials.  Table  3  presents 
the  calculated  results  for  several  materials  which  are  of  interest. 

It  is  seen  from  Table  3,  that  nearly  100#  of  an  incident  ultrasonic  wave 
is  reflected  at  a  gas-solid  interface. 
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Table  3 


Percentages  of  Reflected  Energy  at 

the 

Interface 

of  Two 

Different  Materials 

Medium  I 

Medium 

II 

ri 

Synthetic  Rubber 

86 .711* 

Glass 

20.267 

Steel 

Water 

20.267 

Air 

99-996 

Glass 

68.595 

Synthetic  Rubber 

water 

68.595 

Air 

99.883 

Water 

0 

Glass 

Air 

99.989 

Water 

Air 

99.989 

Attenuation  of  an  Ultrasonic  Wave 

The  overall  attenuation  of  an  ultrasonic  wave  in  a  solid  is  a  complex 
phenomenon  which  depends  on  such  factors  as  bond  breaking,  scattering, 
and  thermal  conduction  of  the  crystals  or  particles  of  the  medium.  Because 
of  the  aforementioned  factor  the  attenuation  phenomenon  results  in  an 
exponential  decay  of  the  wave  intensity  as  it  propagates  in  a  medium. 
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Furthermore,  with  each  propagation  medium  there  is  associated  a 
characteristic  decay  constant,  a.  Thus 


I 

x 


I 

o 


-ax 

e 


(27) 


where 

I  is  the  attenuated  intensity  of  the  ultrasonic  wave  at  distance 
x  from  the  origin,  and 

I  is  the  intensity  of  the  ultrasonic  wve  at  the  origin. 

The  absorption  coefficient,  is  especially  large  when  the  mean 
dimension  of  the  small  crystals  of  the  medium  is  of  the  order  of  a  wave¬ 
length  (16).  For  the  majority  of  solids  and  metals  the  damping  of  sound 
may  be  determined  by  the  following  equation  (l6).  Thus 


a  =  Af  ♦  Bf 


4 


(28) 


where  A  and  B  are  empirical  determined  proportionality  constants.  The 
first  term  in  equation  28  predominates  in  the  majority  of  solids  at  low 
frequencies  and  is  due  to  the  presence  of  elastic  hysteresis,  the  second 
term  accounts  for  the  scattering  of  sound  energy  by  fine  particles.  A 
complete  discussion  of  the  frequency  dependence  of  the  absorption  upon 
the  wavelength  and  crystal  dimensions  may  be  found  in  References  16,  17. 


Ultrasonic  Method  of  Burning  Rate  Measurement 
The  relative  position  of  the  receding  burning  surface  of  a  solid 
propellant  sample  ma'r  be  determined  in  general  by  three  different 
methods:  1.  attenus  .on  measurements,  2.  resonance  measurements,  and 
3.  pulse-echo  measurements.  Further,  the  experimenter  has  the  choice 
of  either  propagating  the  ultrasonic  wave  through  the  combustion  gases  or 
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through  the  solid  propellant  sample.  However,  the  measurement  of  either 
the  intensity  or  the  propagation  of  ultrasonic  waves  in  the  combustion 
gas  will  be  affected  if  acoustic  resonance  is  present  in  the  research 
rocket  motor  (l8).  In  addition,  the  acoustic  waves  produced  by  the 
exhaust  nozzle  will  cause  experimental  errors.  Thus,  if  the  wave  is 
propagated  through  the  solid  propellant  sample  and  not  through  its  com¬ 
bustion  products  the  wave  will  not  be  subject  to  the  acoustic  field  of 
the  research  rocket  motor  (refer  to  Table  3). 

Attenuation  Measurements 

The  thickness,  and  consequently  the  burning  rate,  of  a  solid  pro¬ 
pellant  sample  as  a  function  of  time  may  be  measured  by  measuring  the 
relative  intensity  (wave  magnitude)  of  the  reflected  and  incident  waves 
as  predicted  by  equation  27-  However,  experiments  by  Redwood  (19)  have 
indicated  that  the  exponential  decay  of  an  ultrasonic  wave  is  subject 
to  significant  distortions. 

The  distortion  of  the  exponential  decay  may  originate  from  propel¬ 
lant  surfaces  which  are  not  parallel,  from  interferences  between  modes 
of  propagation,  from  non-uniform  films  which  are  employed  to  couple  the 
transducer  to  the  specimen,  and  from  the  unequal  stressing  of  a  specimen 
due  to  thermal  gradients.  Since  such  effects  may  be  present  in  a  burning 
propellant  sample  the  determination  of  the  position  of  the  burning  pro¬ 
pellant  surface  by  intensity  measurement  would  be  subject  to  significant 


errors . 
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Resonance  Measurements 

The  resonance  method  of  determining  the  position  of  the  burning 
surface  of  a  solid  propellant  sample  Involves  sending  a  continuous 
ultrasonic  vave  through  the  propellant  sample  A  resonance  of  the 
reflected  and  incident  wave  is  obtained  when  the  thickness  of  the  pro¬ 
pellant  sample  is  an  integral  number  of  half  wavelengths.  However,  in 
order  to  obtain  a  large  number  of  resonances  for  a  small  variation  in 
the  sample  thickness  the  wavelength  must  be  small  ( i  e  the  frequency 
must  be  high)  At  high  frequencies  the  deterninat ion  of  the  resonances 
in  a  medium  is  difficult  because  of  harmonic  resonances  and  highly  damped 
resonance  regions  Thus,  the  resonance  method  is  not  an  attractive  method 
for  burning  rate  measurements 

Pulse- Echo  Measurements 

Figure  7  illustrates  schematically  the  pulse-echo  technique  for 
direct  burning  rate  measurement  of  a  solid  propellant  sample  The 
technique  involves  measuring  the  time  required  for  sonic  pulse  to  travel 
from  the  transducer  to  the  burning  surface  and  then  back  (the  echo)  to 
the  transducer  The  input  pulse,  the  echo  pulse,  and  the  time  elapse 
between  the  two  pulses  can  either  be  displayed  by  an  oscilloscope  or 
recorded  by  an  oscillograph 

As  shown  in  Fig.  7.  a  time  delay  spacer  which  has  a  specific 
acoustic  impedence  similar  to  that  of  the  propellant,  is  inserted  between 
the  ultrasonic  transducer  and  the  propellant  sample.  The  function  of  the 
spacer  is  to  increase  the  elapsed  time  between  the  pulse  and  its  echo, 
thereby  raising  the  precision  and  the  accuracy  of  the  measurement. 


') 


a 


FIG.  7  PULSE-  ECHO  MEASUREMENT  TECHNIQUE 
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The  spacer  will  also  protect  the  ultrasonic  transducer  from  the  high 
temperatures  that  occur  as  the  burning  surface  approaches  the  transducer. 

A  thin  film  of  oil  is  to  be  applied  to  the  interface  of  the  ultra¬ 
sonic  transducer  and  the  time  delay  spacer  to  ensure  a  region  of  high 
transml®ivity.  The  propellant  sample  will  be  bonded  to  the  time  delay 
spacer  and  the  propellant  holder. 

Figure  8  presents  the  block  diagram  of  such  a  system.  A  pulse 
from  the  pulse  generator  propagates  to  the  ultrasonic  transducer  and  is 
simultaneously  recorded  by  an  oscillograph  and  displayed  by  an  oscilloscope. 
The  ultrasonic  pulse  wave  travels  through  the  propellant  sample  and  is 
reflected  at  the  burning  surface,  and  the  reflected  wave  then  travels 
back  to  the  transducer,  which  also  functions  as  an  ultrasonic  detector. 

A  part  of  the  energy,  however,  is  again  reflected  at  that  interface  and 
thus  multiple  reflection  occurs.  Consequently,  the  transducer  receives 
not  only  the  desired  reflected  wave  but  also  the  multi-reflected  waves 
and  the  waves  reflected  by  the  sides  of  the  propellant  sample.  But  the 
amplitude  of  the  multi- reflected  waves  will  be  subject  to  attenuation 
in  the  propellant  sample.  By  proper  design,  a  Schmitt  trigger  apparatus 
may  be  employed  for  discriminating  between  the  first  wave  reflection  and 
the  multi- reflected  waves. 

Figure  9  presents  the  output  of  the  ultrasonic  transducer  as 
displayed  by  an  oscilloscope.  The  trigger  level  is  to  chosen  such  that 
only  the  amplitude  of  the  imput  pulse  and  its  first  reflection  are  above 
the  triggering  level  The  'input  pulse  and  its  reflection  are  then 
recorded  by  an  oscillograph.  Since  the  speed  of  the  oscillograph  and  the 
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FIG.  8  BLOCK  DIAGRAM  OF  THE  ULTRASONIC  SYSTEM 
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ULTRASONIC  TRANSDUCER  OUTPUT 


FIG.  9  ULTRASONIC  TRANSDUCER  OUTPUT 
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time  sweep  of  the  oscillograph  are  known,  the  span  between  these  two 
pulses  is  directly  proportional  to  the  travel  time  of  the  pulse.  By 
knowing  the  propagation  velocity  of  the  ultrasonic  wave  in  the  propel¬ 
lant  sample,  the  burning  rate  can  be  calculated  by  the  following 
equation.  Thus 

r  =  f  As  =  f  VY=^fV^r  (29) 

where 

f  =  the  repetition  rate  of  the  generated  pulse  in  cps, 

A  s  -  the  thickness  change  during  the  time  interval  in  inches, 

V  =  the  propagation  velocity  in  in/sec., 

At  =  the  traveling  time,  in  sec., 

A  L  =  the  difference  in  the  pulse-echo  spans  of  two  successive 
input  pulses,  and 

C  =  the  recording  speed  in  in/sec. 

A  Modification  of  the  Ultrasonic  System 
A  possible  modification  of  this  system  would  be  to  use  a  ramp  unit, 
as  shown  in  Fig.  10,  for  measuring  the  traveling  time  of  the  reflected 
wave.  The  ramp  function,  i.c.,  a  voltage  having  a  value  which  increases 
linearly  with  time,  is  started  at  the  same  time  that  the  input  pulse  is 
generated.  The  voltage  increases  linearly  until  it  is  stopped  by  the 
echo  pulse.  For  each  interval  of  time  between  the  input  and  the  echo 
pulse,  the  peak  voltage  will  be  a  direct  measurement  of  the  time  interval, 
and  the  latter  can  be  converted  to  a  voltage  output  with  respect  to  the 
time.  The  time  interval,  which  is  proportional  to  the  burning  rate,  may 
be  recorded  directly  by  an  oscillograph. 


FIG.  10  M00IFIED  ULTRASONIC  SYSTEM 
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Status  of  the  Ultrasonic  Investigation 

Currently,  experiments  are  being  conducted  for  determining  an  optimum 
ultrasonic  frequency  for  non- burning  samples  of  both  aluminized  and  non- 
aluminized  composite  propellants.  An  optimum  ultrasonic  frequency  for  each 
type  of  propellant  will  provide  for  the  desired  attenuation  in  the  ampli¬ 
tude  of  the  echo-pulse. 

After  the  optimum  ultrasonic  frequency  has  been  determined,  small 
samples  of  propellant  will  be  fired  in  a  two-dimensional  research  rocket 
motor  which  will  incorporate  the  ultrasonic  transducer.  The  data  obtained 
by  the  ultrasonic  technique  will  be  compared  with  known  burning  rate  data 
for  the  subject  propellant. 

Experiments  will  be  conducted  to  determine  the  effects,  if  any,  of 
the  interaction  of  the  propellant  combustion  zone  and  the  ultrasonic 
pulses . 
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APPENDIX  A 

NOTATION 


Servomechanism  Technique 

Symbol  Explanation 

E  Feedback  signal  statistical  error 

n  Nominal  counting  rate 

r  Burning  rate 

r  Velocity  of  propellant  positions  shaft 

6 

Root  mean  square  error 

Microwave  Technique 

Symbol  Explanation 

c  Velocity  of  light 

e  Electron  charge 

E  Electric  field  intensity 

f  Microwave  frequency 

f  Mean  frequency  of  quantum  states  1  and  m 
J  m 

h  Plant's  constant 

1  Convective  current 

j  j  - 

k  Boltzmann  constant 


Units 


Count  s/Time 
In/Sec 


Units 

cm/ sec 

coulomb 

Volt/m 

cps 

cps 

erg/ sec 

coulomb/  m‘  sec 
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Symbol 


M 


N 

P 

t 

T 

V 

x 

Z 

Greek 


uT 

P 


Explanation 
Molecular  weight 

Concentration  of  an  ionized  gas  molecule 

Molecular  concentration 

Avogadro's  number 

Gus  pressure 

Time 

Temperature 

Velocity  of  the  electron 


n 


i 


Degree  of  ionization  X  =  — 

o 

Thickness  of  the  combustion  chamber 


a  Absorption  coefficient 

€  Permissivity  of  the  medium 

M  Permeability  of  the  medium 

pQ  Density 

Pe  Charge  density 

T"  Electric  conductivity  of  the  medium 

LfJ  Microwave  radian  frequency 


Plasma  radian  frequency 

_  i  ^  —  V 

Gradient  operator  i  ♦  J  ♦  k  — 


Units 

lb 

1/ cm 

1/ cm^ 

1 / gm-mole 

mm  Hg 

second 

o._ 


ft/j 


cm 


ec 


l/ m 

fa rad/ m 

nenry/m 

lb/ft3 

coulomb/ft3 

mho/ cm 

1/sec 

1/ sec 
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Ultrasonic  Technique 


Symbol 

Explanation 

Units 

C 

Oscillograph  recording  speed 

in/ sec 

f 

Ultrasonic  frequency 

cps 

Ji 

Intensity  of  the  incident  wave 

erg/  cm1' 

Xr 

Intensity  of  the  reflected  wave 

erg/ cm 

I 

o 

Initial  intensity 

erg/ cm2 

I 

X 

Intensity  at  distance  x 

erg/ cm2 

r 

Burning  rate 

in/ sec 

V 

Sonic  velocity 

cm/ see 

X 

Ultrasonic  travel  ing  distance 

cm 

z 

Specific  acoustic  impedence 

gm/cm2sec 

Greek 

a 

Absorption  coefficient 

l/ cm 

b>  s 

Thxness  change  of  the  propellant 

Inch 

bt 

Ultrasonic  traveling  time 

sec 

Difference  in  the  pulse-echo  spans 
of  two  successive  input  pulses 

inch 

inch 

M 

0 

Density  of  the  medium 

g  worn 

APPENDIX  B 


DERIVATION  OF  THF  GENERAL  WAVE  EQUATION 
Maxwell's  equations  for  a  free  space  may  be  expressed  as: 


V-  I  -  Pe/e  (30) 
V*  B  *  0  (31) 
Vx  E  -  -  H  (32) 

VxH=I  +  £  |  (33) 


where  E  is  the  electric  field  intensity,  B  is  the  magnetic  flux  density, 

H  is  the  magnetic  field  intensity,  is  the  change  density,  i  is  the  con¬ 
vective  current,  and  €  is  the  permissivity  of  the  medium. 

A  convenient  vector  identity  is: 

Vx  Vx  E  »  (3U) 

However,  taking  the  curl  of  eqn.  32  yields 

V*  (  Vx  E)  -  Vx  (  -  m|) 

Vx  (  V xE)  !  *  (  V  x  H)  (35) 

Utilizing  eqn.  33,  eqn.  36  may  be  written  as: 

Vx(  VxE)--m£  [i*«  j] 

•  •  M|t  1  •  ^  (36) 

ot 

Thus,  by  eqn.  34  the  general  form  of  the  wove  equation  is: 

V(  V-  E)  -  X?2  E  •  -  vjg  -  f* 

St 


(37) 


